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ABSTRACT 

Aims. The planetary system around the M star Gliese 581 contains at least three close-in potentially 
low-mass planets, Gb 581 c, d, and e. In order to address the question of the habitability of GL 581 d, 
we performed detailed atmospheric modeling studies for several planetary scenarios. 
Methods. A ID radiative-convective model was used to calculate temperature and pressure profiles of 
model atmospheres, assumed to be composed of molecular nitrogen, water, and carbon dioxide. The 
model allows for changing surface pressures caused by evaporation/condensation of water and carbon 
dioxide. Furthermore, the treatment of the energy transport has been improved in the model to account 
in particular for high CO2, high-pressure Super-Earth conditions. 

Results. For four high-pressure scenarios of our study, the resulting surface temperatures were above 
273 K, indicating a potential habitability of the planet. These scenarios include three C02-dominated 
atmospheres (95% CO2 concentration with 5, 10, and 20 bar surface pressure) and a high-pressure 
C02-enriched atmosphere (5% CO2 concentration with 20 bar surface pressure). For all other consid- 
ered scenarios, the calculated GL 581 d surface temperatures were below the freezing point of water, 
suggesting that GL 581 d would not be habitable then. 

The results for our C02-dominated scenarios confirm very recent model results by Wordsworth et al. 
(2010). However, our model calculations imply that also atmospheres that are not C02-dominated (i.e., 
5% vmr instead of 95% vmr) could result in habitable conditions for GL 581 d. 

Key words. Astrobiology, Planets and satellites: atmospheres, Stars: planetary systems, Stars: individual: 
Gliese 581, Planets and satellites: individual: Gliese 581 d, 
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1. Introduction 
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Since the discovery of the first extrasolar plan- 
ets around pulsars ( Wolszczan Sz Frail 19921 ) and 
-sequence stars (jMavor fc Quelozlll995l " more 



mam- 



than 400 planets and planet candidates have been 
detected, most of which are Jupiter-mass objects. 
However, in recent years, detections of smaller plan- 
ets have been announced, and today 21 planets with 
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(minimum) masses below 10 Earth masses (mm) ar e 
known (e.g. lRivera et al.ll2005l . lHoward et al.ll2 009) 
Two of these planets, CoRoT-7 b (ILeeer et al 



20091 ) and GL 1214 b dCharbonneau et ahl l2009T l 



have been discovered by the transit method and 
produce a measurable radial velocity variation, 
hence their true mass and mean density could be 
derived. They are the first two examples of so-called 
"Super-Earths" (1 < m < 10 m®). 

Modeling studies of the atmospheres of 



ter- 



restrial (exo)planets (e.g. 


Kastine et al. 


1993, 


Joshi et al. 1997. Foreet & Pierrehumbert 


1997 


. Pierrehumbert <fe Erlick 1998. Selsis et al. 


2002 


. Seerura et al. 2003. 


Seeura et al. 


2005, 


Grenfell et al.l 


2007bl. Grenfell et al. 


2007al. 


Kitzmann et al. 


2010) have 


aimed at charac- 



terizing the response of the atmospheric system 
to changes in planetary (e.g. atmospheric compo- 
sition, presence of clouds) and stellar parameters 
(e.g. central star type, orbital distance). These 
studies mostly assessed the potential surface habit- 
ability of such planets and predicted their spectral 
appearance. Surface habit ability is commonly re- 
lated to the presence of liquid water on the surface, 
hence temperatures above 273 K, because liquid 
water seems to be the fundamental requirement for 
life as we know it on Earth. 

In view of these modeling activities , the system 



Gliese 581 (GL 581) with four planets (IBqnfils et al. 



20051 . lUdrv et al.ll2007l . lMavor et alJl2009 ) is partic- 
ularly interesting. It hosts at least three potentially 
low-mass, therefore possibly terrestrial planets. The 
one closest to the central star , GL 581 e, was an- 
nounced by Mayor et al. ( 20091 ) and has a minimum 
mass of 1.94 mg with an orbital distance of 0.03 
AU to the star. The two outer low-m ass planets 



GL 5 8 1 c and d, were di scovered by lUdrv et al 
(120071 ). iMavor et all d2009l ) refined the orbital dis- 
tances and minimum masses of these two planets, 
obtaining 5.36 tue (GL 581 c at a distance of 0.07 
AU to the star) and 7.09 m E (GL 581 d at a dis- 
tance of 0.22 AU to the star), respectively. 



Taking 



into 



account 



(LoDez-Morales et al. 




2006, 


Mavor et al. 


2009) 


and dynamical ( 


Beust et al. 2008. 


Mavor et al. 



system most probably lies between 40° and 85°. 
This implies masses no higher than 1.56 times the 



minimum masses, so all three low-mass planets are 
likely to be Super-Earths. 

The habita bility of GL 581 c and d has been inves- 



tigated by ISelsis et al.1 (|2007l ) and Ivon Bloh et al 
(12007a!). ba s ed on the discovery da ta published by 



Udrv et al. (|2007l ). ^n their study, Ivon Bloh et al 



(|2007al ) applied a box-model of the Earth, which 
incorporates the carbonate-silicate cycle, requiring 
a surface reservoir of liquid water and a tectoni- 
cally active planet. It includes the exchange of CO2 
between the mantle and crust of the planet and its 
atmosphere by assuming parameterizations for con- 
tinental growth and spreading rate. The CO2 par- 
tial pressure and the stellar luminosity were related 
to the surface temperature through a simple energy 
balance equation b etween stellar and emit ted ther- 

Bloh et al.ir2007bl for 



mal fluxes (see also 



von 



more 



details) . ISelsis et al.1 (|2007l ). on the other hand, used 
previous model results for the habitable zone (H Z) 
from the seminal study of iKasting et al. (Il993l ) to 
estimate the boundarie s of the HZ in the G L 581 
system. The results of Kasting et al. ( 19931 ) were 
obtained with a ID radiative-convective model. 



Furthermore. ISelsis et al.1 (120071) used mod el results 
for early Mars from Mischna et al. ( 2000l ) to illus- 
trate the uncertainty in the outer limit of the HZ 
becaus e of the (pos si ble) p resence o f CO2 clouds. 



Both Selsis et al. 



von Bloh et al 



(|2007l ) and 

()2007al ) concluded that the inner planet GL 581 
c is unlikely to be habitable, because it is closer 
to the star than the inner boundary of the HZ, 
whereas the outer planet, GL 581 d, ju st might be 



habitable. Based on t he ca lculations bv ISelsis et al 
()2007l ). lMavor et al.l (120091 ) concluded that GL 581 
d lies in the HZ of its central star, considering that 
the refined orbit means it receives more than 30 % 
more stellar e nergy than previous l y tho ught. 
Very recently, Wordsworth et al. ( 2010l ) presented 
ID modeling studies of high CO2 atmospheres of 
GL 581 d, varying the CO2 pressure and other 
parameters such as relative humidity and surface 
albedo. They found habitable surface conditions 
with CO2 parti al pressures of >5 bar. I n a simi- 
lar approach to IWordsworth et al. (|2O10l ). we ap- 
plied a different ID radiative-convective model to 
the atmosphere of GL 581 d considering different 
appropriate planetary scenarios. By calculating at- 
mospheric pressure, temperature, and water pro- 
files, a reasonable range of planetary surface condi- 
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tions has been investigated with regard to potential 
habitability. 

The paper is organized as follows: Sect. [2] describes 
the properties of the planetary system GL 581. The 
model used is described in Sect. El A description of 
the model input parameters and of the runs per- 
formed is given in Sect. [H Results are described 
and discussed in Sect.[5j We give our conclusions in 
Sect. El 



D -18_ 




IUE 
data 



100 



wavelength [ftrn 



2. The planetary system GL 581 



2.1. Main properties of the star GL 581 



GL 581 is a quiet M3 star (jlldrv et al.ll2007h . Its 
main parameters are summarized in Table [TJ 



Table 1. Main properties of GL 581 



Property 



Value 



Reference 



Type M3 Udrv et al. f 20071 

Mass 0.31 M Bonfils et al. ('2005') 

Luminosity 0.013 L Udrv et al. f 20071 . Bonfils et al. (20051 

0.0135 L Selsis et al. f 20071 ~ 

Radius 0.38 R Lacv (19771 . Johnson fc Wright 1 19831 

0.29 R Chabrier k Baraffe (20001 

T e fl 3190 K calculated with the Stefan-Boltzmann law 

3249 K Johnson h Wright (19831 

3760 K Butler et al. (20061 

[Fe/H] -0.10 Johnson fc Adds f 20091 

-0.25 Bonfils et al. (20051 

distance 6.27 pc Butler et al. (2006 1 

6.53 pc Johnson fc Wright (19831 . Lacv (19771 



In order to perform proper atmospheric modeling 
studies, we need a spectrum of GL 581. Therefore, 
we derived a high-resolution spectrum (Fig. [1]) from 
two sources, i.e. an UV spectrum measured by the 
IUE (International Ultraviolet Explorer) satellite 
i n 1987 and a syntheti c Nextgen model spectrum 
dHauschildt et alJll99d ). 



Fig. 1. Derived high-resolution spectrum of GL 581 
(data sources as indicated). 

The model parameters of the available NextGen 
model spectrum used here are T e g=3200 K, log 
g= 4.5, [Fe/H]=0.0 and i? G L58i=0.38 R & . The 
synthetic spectrum was calculated for wavelengths 
from 0.334-971 /um. It was merged with the IUE 
spectrum from 0.185-0.334 /im. 
Based on the 6 .27 pc distance of GL 581 
(jButler et al.ll200fih . the spectrum was scaled to the 



orbit al distance of GL 581 d (0.22 AU, iMavor et al 



2009h and binned to the spectral intervals required 
for the model code (see Sect. [3j). 
The orbit of GL 581 d could be highly eccentric 
(Table [2]). The mean stellar flux F received over an 
eccentric orbit is given by 



F(a) 



(1 



o2\0.5 



where F(a) is the flux received by the planet at 
th e distance of its semi-maj or axis a. As shown 
bv I Williams fe Pollard! foooi ) with a 3D model ap- 
proach, the overall climate of a planet on an ec- 
centric orbit behaves almost as if the planet were 
constantly receiving this averaged flux F. For GL 
581 d, the application of Eq. ([1]) results in an in- 
crease of the stellar flux of about 8% compared to 
the circular case. 

Assuming that the atmospheric response timescales 
to changes in stellar flux (which are on the order of 
weeks/months) are comparable to the orbital pe- 
riod of the planet (about teo months), F was used 
for the model calculations. Note, though, that this 
approach depends on the properties of the central 
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star and the planetary system. For larger orbital 
distances, for instance, flux variations caused by or- 
bital eccentricities should be taken into account. 



2.2. Properties of the planet GL 581 d 

The semi-major axis is known from the measured 
orbital periods and the mass of the ce ntral star via 



-Asiir f=0-13 (Earth mean, see iKitzmann et al.l 12010 
and Rossow fc Schiffer 19991 ). 

Table [2] summarizes the assumed planetary param- 
eters of GL 581 d. 



Table 2. Planetary parameters of GL 581 d. 



Keple r's third law. This we took from lMavor et al 
<|200gh to be a GL 58id = 0.22 AU. The best fit of the 



Property 



Value 



Reference 



radial velocity data was obtai ned for a highly ec - 
centric or bit with e d =0.38 (see lrvlavor et alJl2009h . 



In order to be cons istent with the studi es of 
Selsis et al.l ((20071 ) and lvon Bloh et all (j2007ah . the 
stated minimum mass was assumed to be the true 
mass m of the planet. The planetary radius r is 
then taken from a theoretica l mass-radius relat ion- 
ship for terrestrial planets (jSotin et al.l 120071 ) . fi- 



7.09 m e Mayor et al. (2QQgj 

0.22 AU Mayor et al. (2009) 

0.38 Mayor et al. (2009*1 

1.71 r ffl after Sotin et al. (2007) 

23.76 ms~ 2 calculated 

0.13 Earth mean, Rossow k Schiffer (1999) 



nally yielding the surface gravity. 

Usually, planets in the HZ of M stars (or closer 
to the star) are assumed to be tidally locked 
in a synchronous rotation, i.e. they rotate with 
the same rate uj as their orbital per iod P (e.g. 



Lammer et al.l 120071 . IScalo et al.1 120071). However 
as already note d by Goldreich Sz Peak _J 19681 ) 



Correia et al 



Mass 

Orbital distance 
Eccentricity 
Radius 
Gravity 
Surface albedo 



3. Atmospheric model 

3.1. General model description 

A ID, cloud-free radiative-convective column model 
was used for the calculation of the atmospheric 
structure, i.e. the temperature and pressure profiles 
for different GL 581 d scenarios. 
The mod el is originally based on the one de- 
scribed by Kasting et al. ( 1984al ) and Kasting et al 



( 20081 ). planets that show significant 
eccentricities - possibly like GL 581 d - are unlikely 
to become locked in a l:l-resonance. Therefore, GL 
581 d is not considered to rotate synchronously in 
this study. 

An important model parameter for ID radiative- 
convective models is the surface albedo. One ap- 
proach used in these models i s to adjust the sur- 
face albedo in the mo d el (e.g.. | Kasting et al.lll993l . 



e.g. 



GrieBmeier et al.| ([20051) , | Levrard et al. | ([20071) and di984blh Fu r ther d evelopments ar e described by 



Kastinel (119881). iKastind (Il99ll ). iKasting et al 



1993), iMischna et al.1 (|2000l ) and IPavlov et al 



(|2000i ). Additional updates of the thermal radia- 
tio n scheme of t he m odel have been introduced 
by Segura et al. ( 20031 ). The model version used 



here is based on the version of von Paris et al 



Grenfell et al. 2007b. Goldblatt et al 



2009 



so 

that the models reproduce prescribed reference sce- 
narios (e.g. modern Earth with a surface tempera- 
ture of 288 K). Therefore, the values of the model 
surface albedos largely depend on the impact of 
clouds on the surface temperature. This impact, 
however, is unknown for Super-Earths. Hence, in- 
stead of an arbitrary choice of the surface albedo 
used in the tuning of ID models, we adopted the 
measured surface albedo of Earth here. Thus, the 
unconstrained effect of clouds in Super-Earth at- 
mospheres is explicitly excluded from our cloud-free 
simulations. The surface albedo value used here is 



(|2008l ) where a more detailed model description is 
given. The current model uses the radiative trans- 
fer scheme MRA C in the IR (see the introduction 
of this scheme in von Paris et al. 20081 ). The water 
profile in the model is calcul ated based on the rela- 
tive h umidity distribution of lManabe &: Wetherald 
(|1967h . Above the cold trap, the water profile is set 
to an isoprofile at the cold trap value. 
The model considers N2, H2O, and CO2 as atmo- 
spheric gases. Other radiative trace species might 
be present in the atmospheres of exoplanets (e.g. 
SO2, O3, or CH4), which could alter the radiative 
budget significantly. But because the presence of 
these gases highly depends on the planetary sce- 
nario (e.g., outgassing, volcanism, formation, bio- 
sphere), our model atmospheres are restricted to 
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the two most important greenhouse gases of the 
Earth's atmosphere (H2O and CO2), using N2 as 
additional background gas. This is based on the 
observation that in the solar system all terrestrial 
atmospheres (Venus, Earth, Mars, Titan) contain 
significant amounts of N2 . 

Temperature profiles from the surface up to the 
mid-mesosphere are calculated by solving the equa- 
tion of radiative transfer and performing convective 
adjustment, if necessary. The convective lapse rate 
is assumed to be adiabatic. Water profiles are cal- 
cula ted by assuming a fixed rela tive humidity pro- 
file dManabe fc Wetheraldl fl967h through the tro- 
posphere. 

As a globally averaged ID column model, the model 
does not incorporate dynamical processes such as 
winds and latitudinal energy transport. Thus po- ^" 
tentially important 3D effects are not included. 
For example, the effect of atmospheric collapse and 
global energy re-distri bution on tidally lo cked plan- 
ets wa s investigated b ylJoshi et al.1 ( 19971) a nd lJoshi 
(|2003h . The study of ISpieeel et al.1 (|200Sl ) showed 
that slowly rotating planets feature a very effi- 
cient energy redistribution because of the weaker 
Coriolis force (compared to Earth). Hence a ID 
model could give a reasonably accurate picture of 
global mean temperatures (and therefore, habitabil- 
ity) on planets orbiting M stars. In addition to ro- 
tation, also atmospheric optical depth and mass of 
a planet are important for assessing temperature 
gradients (jJoshi et al,lll997l ). Therefore, 3D atmo- 
spheric modeling of GL 581 d is needed in the fu- 
ture to investigate its potential habitability in more 3. 
detail. 

3.2. Computational details 

In this work, we employed an impr oved version of ^" 
the model of von Paris et al. ( 20081 ). 



1. In the present study, CO2 condensation is in- 
cluded for the calculation of the adiabatic lapse 
rate in the c onvective regim e acc ording to the 



treatm ent by iKastind (|199ll ) and iKasting et al. 



( 19931 ). Th e saturation vap or pressure of CO2 is 
taken from I Ambrose! ( 19561 ) . CO2 condensation is 
generally assumed to occur when the atmosphere 
is (super-)saturated with respect to CO2. This 



criterion is expressed by the super saturation ra- 
tio Sgat, which should be higher than unity: 



PCO2 



P. 



Ssat > 1 



sat,C0 2 



where pco 2 ls t ne partial CO2 pressure and 
-Psat,co 2 the saturation vapor pressure of CO2. 
The numerical value of 5 sat used here is 
taken to be 1 - 34, b ased on measurements by 
Glandorf et all l|2002h . This implies that the at- 
mosphere has super-saturated with respect to 
CO2 for CO2 c onden sation to occur . Note that 
Kastind (|199ll ) and IKasting et all (|1993l ) as- 
sumed S , sat =l, i.e. the model atmosphere was as- 
sumed to be saturated. 

The calculations now also include Rayleigh scat- 
tering by water vapor (cr myi H 2 0] m [cm~ 2 ]) cal- 
culated via 



<Tray,H 2 o(A) = 4.577- 10 



-21 



6 + 3-D 



6-7-DJ A 4 



7! (3) 



where D is the depolarization ratio, r the refrac- 
tivity and A the wavelength in 11 m. The nume rical 
factor 4.577- 10~ 21 is taken from [Men] ()l973h . 

This work assu mes D = 0.17 from 

Marshall fc Smith! (jl990l ). The refractivit y 
is calculated as r = 0.85 • rd rya ir ( Edlen 19661 ). 



The refractivity of dry air is obtained from an 
approximate formula given bv lBucholtzl (Il995h . 



The heat capacity of water vapor has been in- 
cluded in the calcu lations based on the so- called 
Shomate equatio n (Parks Sz Shomate 1940l ) with 
parameters from IChasd (| 19981 ). 



One critical issue for the IR radiative trans- 
fer scheme used here is the so-called bi- 
nary species parameter 77 (Mlawer et al. 19971 . 



Colanrete fc Toon! 12003 . Ivon Paris et all 120081 ). 



In previous studies, the variation of absorption 
cross sections of different mixtures was assumed 
to be linear according to the relative concen- 
trations of the major absorbers (here, H2O and 
C02). For atmospheres where the relative con- 
centrations of water and carbon dioxide vary by 
several orders of magnitude, this approach de- 
pends on the assumed reference ratio. Hence, a 
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logarithmic 7] grid has been introduced where rel- 
ative concentrations are allowed to vary by 15 
orders of magnitude. 

5. The IR H2O contin uum calculation used in 



profiles. Further comparative studies with the im- 
proved model are presented in Appendix 1X1 



von Paris et al. ( 20081 ) has b een replaced by the 



CKD continuum formulat ion (jClough et al.lll989l . 
Schreier &: Bottger 120031 ). In this way, the 



con- 
tinuum is no longer restricted to the atmospheric 
window region (between 8 and 12 /im) and in- 
corporates contributions from H2O self as well 
as foreign continua and, for the first time in the 
context of exoplanet ID models, the CO2 foreign 
continuum in the IR. 

6. The temperature-pressure grid for the absorption 
coefficient calculations in the IR has been ex- 
tended to T=100 K and p=1000 bar This was 
done as described by Ivon Paris et al.l (120081'). u s- 
ing line parameters from Rothman et al. ( 19951 ). 



7. The model has been further modified to account 
for changing pressure levels during the calcula- 
tions. Under some conditions described here, the 
atmospheric content of water or carbon dioxide is 
controlled by the surface temperature caused by 
evaporation or condensation. The surface pres- 
sure is thus determined by the surface tempera- 
ture and must be re-calculated after each model 
time step to adjust to the newly calculated value 
of the surface temperature. The model pressure 
grid is therefore no longer fixed, but dynamically 
determined between two consecutive iterations in 
the model. 

The IR radiative transfer updates mentioned above 
were tested against high-resolution line-by-line (lbl) 
calculation s with the MIRART ra diative trans- 
fer code ( Schreier fc Bottger 2003), which was 
found to compare well with other lbl codes (e.g., 
Melsheimer et aI.ll2Q05h . 

The total IR flux was within 5% of the lbl values 
for all atmospheric levels. Fluxes for each spectral 
interval of the radiative transfer code were usu- 
ally within 10% of the lbl values, except for a few 
near-IR bands in the middle to upper atmosphere. 
However, these bands do not contribute strongly 
to the cooling rate (less than 10" 6 ) or the strato- 
spheric energy budget, hence these deviations are 
insignificant for the calculation of the temperature 



4. Model input for GL 581 and performed 
runs 

Important properties characterizing planetary at- 
mospheres are the total surface pressure and the 
composition of the atmosphere. These determine 
the greenhouse effect and the volatile reservoir, 
hence the habitability of a planet. Atmospheric 
pressure and composition are determined by the ac- 
cretion and outgassing history, whether from the 
interior or from impacts (e.g. Pepin 199ll ). at- 
mospheric loss because of escape to s pace (e.g. 
Kulikov et all 120071 . iKulikov et all bood ) or incor- 



poration into planetary surface material. But these 
processes are not known for the GL 581 system. 
Hence, in order to study the influence of atmo- 
spheric properties on the potential habitability of 
GL 581 d, we performed a parameter study consid- 
ering different atmospheric scenarios. 

We assumed that the initial water delivered to the 
planet via outgassing and impacts has been re- 
tained, hence a reservoir of liquid water is present 
on the planet ary surface . This a ssumption was 



also m ade b y ISelsis et al. (l2007h. Ivon Bloh et al 



(|2007ah and IWordsworth et al.l <|2O10h for the GL 



581 planetary system. 

The total surface pressure (1, 2, 5, 10, 20 bar) and 
CO2 concentration (0.95, 0.05 and 355 ppm vmr, 
respectively) were varied. N2 is assumed to act as 
a filling background gas. The range of surface pres- 
sures was chosen to represent scenarios adopted in 
the literature for early Earth and early Mars in 
terms of atmospheric column density. The values of 
CO2 concentrations are chosen to represent modern 
Earth (355 ppm CO2 concentration), early Earth 
(5% C0 2 ) and Mars or Venus (95% C0 2 ) condi- 
tions, i.e. representative values of terrestrial atmo- 
spheres within the solar system. Table [3] summa- 
rizes the set of runs for GL 581 d performed here. 
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Fig. 2. Temperature-pressure profiles for the low CO2 (355 ppm CO2) runs of Table [3l The equilibrium 
temperature of the planet (dashed) and melting temperature of water (dotted) are indicated as vertical 
lines. 



Table 3. Atmospheric scenarios for GL 581 d 



Set p [bar] CO2 vmr 

Gl (low C0 2 ) 1,2,5,10,20 3.55 • 10" 4 

G2 (medium C0 2 ) 1,2,5,10,20 0.05 

G3 (high C0 2 ) 1,2,5,10,20 0.95 



5. Results and discussion 

5.1. Temperature profiles 

Figure[2]shows the temperature-pressure profiles for 
the low C0 2 set of runs (Gl) of Table El For all 
runs the surface temperatures are far lower than 
the freezing point of water (273 K, indicated by the 
dotted vertical line in Fig. |2J), which implies that 
liquid water is not present on the surface of the 
planet. Upon increasing the surface pressure from 
1 bar to 20 bar, the surface temperatures increase 
by 4.2 K (from 203.6 K to 208.9 K). 
Also shown in Fig. [2]is the equilibrium temperature 
T eq of the planet. It is about 197 K, hence habitable 
conditions would require a greenhouse effect (GHE) 
of 76 K. As can be inferred from Fig. [21 the actual 
GHE in the model atmospheres is about 6 to 12 



K, which is much less than the approximately 30 K 
GHE provided by the modern Earth atmosphere. 
An interesting feature of the temperature profiles 
of the low CO2 case is the absence of a convective 
troposphere for all Gl runs. These atmospheres are 
all in radiative equilibrium. 

Figure [3] shows the temperature-pressure profiles for 
the medium C0 2 set of runs (G2) of Table [3l For 
the 20 bar run, the value of the surface temperature 
was 313.3 K, hence exceeded the freezing point of 
water. Thus, in this particular scenario GL 581 d 
might be habitable. The other runs with 1, 2, 5, 
and 10 bar though all showed surface temperatures 
below 273 K. Note, however, that for the 10 bar run 
the obtained surface temperature is 268.1 K, which 
is fairly close to the freezing point of water. 
In contrast to the low CO2 case, the increase of sur- 
face pressure from 1 to 20 bar has a huge effect on 
surface temperature, which increased by about 105 
K. This is caused by an enhanced absorption of stel- 
lar radiation by C0 2 and H 2 near-IR absorption 
bands and a massive greenhouse effect, which re- 
sults in the onset of convection for runs with higher 
surface pressures. 

Figure 0] illustrates this effect, showing net (i.e., 
frequency-integrated) stellar and thermal down- 
ward (Ffi) and upward (F u ) fluxes for the medium 
CO2 (G2) 20 bar run. Much of the incoming stellar 
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Fig. 3. Temperature-pressure profiles for the medium CO2 (5 % CO2) runs of Tabled The equilibrium 
temperature of the planet (dashed) and freezing temperature of water (dotted) are indicated as vertical 
lines. 

Stellar fluxes Thermal fluxes 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 j l 1 1 1 1 1 1 ; 




Flux [Wm~ 2 ] Flux [Wm~ 2 ] 



Fig. 4. Net fluxes of the 20 bar run of set G2 (medium CO2): Stellar (left) and thermal (right). 
Downwelling (solid) and upwelling (dotted) fluxes are shown. 



radiation (~ 70 %) is absorbed by CO2 and wa- 
ter in the lower atmosphere, as illustrated by the 
left panel in Fig. 01 The difference C r = F u — 
for the thermal fluxes is the radiative cooling. A 
strong GHE is indicated by a low value of C r . In 
the lower atmosphere of the 20 bar run, both ther- 
mal components balance each other, as can be seen 
in Fig. HI meaning that the thermal radiation is effi- 



ciently trapped in the atmosphere. The value of C r 
in the bottom atmosphere layer is about 5 Wm -2 , 
which corresponds to about 1 % of the surface emis- 
sion (ub T^ ur{ 540 Wm" 2 , <jb Stefan's constant, 
^surf=313 K surface temperature). On Earth, the 
value of C r is about 70 Wm -2 , about 20 % of the 
surface emission. 
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The atmospheric structure of the medium CO2 runs 
in Fig. is very different from the low CO2 runs 
in Fig. [2j Firstly, the temperature inversion in the 
upper atmosphere of the G2 runs at pressures below 
10 mbar is much more pronounced (~ 30 K) than 
for the low CO2 runs. This is owing to the strong 
absorption of stellar radiation by CO2 and H2O in 
the near-IR bands (at 2, 2.7 and 4.3 fim). 
Secondly, convective tropospheres develop for the 
5, 10, and 20 bar runs. These extend to an altitude 
of about 0.5-1 surface scale heights (about 2-5 bar) 
above the surface, which is comparable to the tropo- 
sphere extension on present Earth. Hence, the lapse 
rate is much steeper in the medium CO2 cases than 
in the low CO2 cases. For example, the convective 
lapse rate in the 20 bar medium CO2 case is about 
22 K km -1 near the surface (the dry adiabatic lapse 
rate varies linearly with gravity) , whereas the radia- 
tive lapse rate for the 20 bar low CO2 run is only 
about 9 K km -1 . The value of 22 K km -1 is very 
close to the dry adiabatic lapse rate of 23 K km -1 
in the medium CO2 case. Despite the high surface 
temperature of 313 K and a corresponding partial 
pressure of water of about 70 mbar, the water con- 
centrations are only of the order on 10~ 3 vmr near 
the surface. Hence, the lapse rate is close to the 
dry adiabat (see discussion in Ingersoll 19691 ). even 
if appreciable amounts of water are present in the 
atmosphere. 

In Fig. [5j the temperature-pressure profiles for the 
high CO2 set of runs (G3) are shown. Except for 
the low-pressure runs with 1 and 2 bar surface 
pressure, all scenarios showed surface temperatures 
above 273 K. Hence, the results suggest that CO2- 
rich atmospheres may imply habitable conditions 
on GL 581 d. 

The atmospheric structure in the high-C02 case is 
different from those in the low and medium CO2 
cases. Even the 1 and 2 bar runs now show convec- 
tive tropospheres, albeit not very extended ones. 
More massive tropospheres develop in the 5 bar, 10 
bar, and 20 bar runs owing to the onset of CO2 con- 
densation in the middle atmosphere, as indicated 
by the CO2 condensation curve in Fig. The re- 
spective tropopauses are located 3-5 surface scale 
heights above the planetary surface. On present 
Earth, this would correspond to tropopause heights 
of about 20-40 km compared to the roughly 10 km 
tropopause altitude observed today. 



The tropospheres in the 5, 10, and 20 bar high CO2 
runs are divided into two regimes, an upper tropo- 
sphere with CO2 condensation (see CO2 condensa- 
tion curve in Fig. [SJ and a lower troposphere with 
H2O condensation. This temperature structure is 
comparable to atmospheric structures c alculated 
for models of the early Mars atmosphere ( Kasting] 
199ll . iMischna et allbood . IColaprete fc Toonll2003h 



From Fig.[5]it is also evident that temperature pro- 
files are similar to each other for the 5, 10, and 
20 bar runs. The temperature profiles of these runs 
differ by about 2-5 K at equal pressures, as is the 
case for the 10 and 20 bar medium CO2 runs (Fig. 
[3j). This is because of two reasons. Firstly, the lower 
atmospheres of the 5, 10, and 20 bar high-C02 runs 
become optically thick for thermal radiation, as is 
observed for the medium CO2 10 and 20 bar runs. 
Secondly, an increase in surface pressure does not 
result in an increase in planetary albedo because 
GL 581 emits only a negligible amount of radiation 
in the spectral range where the Rayleigh cross sec- 
tion is largest. Therefore, the global energy balance 
is no longer affected by the lower part of the model 
atmosphere, hence temperature profiles do not dif- 
fer significantly. 

A widely used assumption for the assessment of 
the HZ is t he so-called m a ximu m greenhouse, in- 
troduced by Kasting et al. ( 19931 ) for an Earth- like 
planet around the Sun. With increasing amounts 
of CO2 in the atmosphere of a planet located near 
the outer boundary of the HZ, the GHE will be- 
come more and more saturated, i.e. the optical 
depth is near or higher than unity for all CO2 
bands. Then, any further increase of CO2 will 
only increase Rayleigh scattering, hence increase 
the planetary albedo. Thus, for increasing CO2 
partial pressure, surface temperatures will show a 
maximum and then decrease with increasing CO2 
pressu re. The same behavi or was found for early 
Mars (jMischna et al.1 12000| ). The different behav- 



ior fbr_oj : u^Jin2hCO^runs co mpared to the result s 
of iKasting et all (|l993T ) and IMischna et all f|2000l ) 
is illustrated in Fig. [6) They are summarized in 
terms of atmospheric column density and normal- 
ized surface temperature T sur f/T cq where T sur f is 
the surface temperature of the planet. As is obvi- 
ous from this figure, the high-C02 runs from this 
work do not feature a maximum greenhouse effect. 
However, a maximum greenhouse effect will proba- 
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Fig. 5. Temperature-pressure profiles for set G3 of Table [3] (95% CO2). The CO2 condensation curve 
(with 5 sa t=l-34 in Eq. [2|) is shown in light gray. 
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Fig. 6. Maximum greenhouse effect: Response of surface temperature to increase of atmospheric CO2 
content. 



bly occur at still higher pressures, hence column 
densities (not investigated here). The absence of 
a maximum greenhouse effect for the runs of this 
work is caused by four factors. 

Firstly, GL 581 emits much more radiation in the 
near- to mid-IR, and less in the visible than the Sun. 
Consequently, the contribution of Rayleigh scatter- 
ing to the planetary albedo is much less efficient for 
planets around GL 581 than around the Sun be- 



cause of the A _4 -dependence of the Rayleigh scat- 
tering cross section. Secondly, the stronger near- 
IR emission of GL 581 leads to more heating 
by near-IR absorption ban ds of H2O and CO q - 



Thirdly, the simulations of iKasting et al.l ()1993h 



were done at constant surface temperatures of 273 
K (hence, constant partial pressure of 6.5 mbar wa- 
ter), which neglects the positive feedback provided 
by increased water vapor at higher surface tempera- 
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tures. Fourthly, because of the higher gravity of GL 
581 d compared to Earth, the same column amount 
of CO2 (i.e., x coordinate in Fig. [6]) is reached at 
much higher pressures, e.g. 20 bar on GL 581 d com- 
pared to 8 bar on Earth. The pressure broadening of 
absorption lines then leads to an enhanced absorp- 
tion in the line wings. Because the line centers are 
usually optically thick, a higher absorption coeffi- 
cient in the line wings can significantly increase the 
overall absorption of radiation. It was recently sug- 
gested that this kind of behavior could have helped 
to warm the early Earth b y invoking higher N2 par - 
tial pressures than today (jGoldblatt et al.l l2009al ) . 
It is also mainly responsible for the fact that surface 
temperatures for the 1 bar high CO2 case are about 
85 K lower than for the 20 bar medium CO2 case, al- 
though both atmospheres contain the same amount 
of C02- When interpreting this result, one has to 
bear in mind the potential uncertainties related to 
the radiative transfer in high-pressure atmospheres 
(see sensitivity study in Appendix iBj). 



5.2. Comparison with other studies of GL 581 d 



For several model scenarios (5, 10, and 20 bar high 
CO2 and 20 bar medium CO2 runs) of this work, 
surface temperatures were found to be above 273 
K, i.e. these results imply habitable surface condi- 
tions on GL 581 d. In all other scenarios, GL 581 
d was found to be uninhabitable with surface tem- 
peratures below 273 K. 



The results are summarized in Fig. [7] in the con- 
sidered parameter space, i.e. the surface pressure - 
CO2 concentration plane. 



Yet with regard to the uncertainties in the radiative 
transfer associated with CO2 continuum absorption 
(see the sensitivity study for surface temperature 
in Appendix [B]), the modeling results regarding the 
habitability must of course be treated with caution. 
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100.0 



Fig. 7. Assumed model atmospheres (Table [3]) for 
GL 581 d: Habitable (filled circles) and uninhabit- 
able (open circles) scenarios. 



Still, our results imply that, given plausible Venus- 
or early Earth-like atmospheric compositions and 
taking into account reasonable uncertainties in the 
radiative transfer formulations, GL 581 d can be 
classified as a habitable planet. 

(|2007ah concluded 



The study of Ivon Bloh et al 



that GL 581 d represents a habitable planet even 
for relatively low CO2 partial pressures of about 
4-5 bar. Taking into account the more than 30 % 
increase in insolation for GL 5 81 d owing to the re - 
vision of orbit al parameters by Mayor et al. ( 20091 ). 



the results of Ivon Bloh et al.l (|2007al ) imply that 
GL 581 d could indeed be habitable for even less 
dense CO2 atmospheres. The findi ngs of our study 



do not support these implications of lvon Bloh et al 
(|2007ah though. 



This disagreement is partly owed to the parame- 
terizations and empirical criteria employed. For ex- 
ample, the parametrization of the relation between 
CO2 partial pressure, pla netary albedo, and sur - 

" d2007ah 
dl9£ " 



face temperature use d bv Ivon Bloh e t al.l 
is based on work by Williams Sz Kastind (j 19971 ). 



Budvkol (jl995h and iGhamberlainl (|l98ft ). These 
studies were done for Earth, orbiting around the 
Sun, and thus did not account for different central 
stars with a different spectral energy distribution 
or more massive terrestrial planets. In contrast, the 
radiative-convective model used here takes these ef- 
fects into account. 
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The study of ISelsis et al.l (|2007l ) concluded that 
a pure CO2 atmosphere without additional green- 
house gases is unlikely to provide habitable con- 
ditions on GL 581 d. Like in the study of 
von Bloh et al.1 l|2007al ). they did this for the then 



known orbital parameters of GL 581 d. ISelsis et al 



also provide a discussion of the uncertainties 
of the limits of the outer HZ in view of early Mars. 
From this discussion, GL 581 d could still be in 
the HZ when CO2 clouds or additional greenhouse 
gases are taken into account . Note that, as already 
stated by ISelsis et al.1 (|2007l ) , the use of parameter- 
iza tions of the ou t er bo undary of the HZ provided 
by Kasting et al. ( 19931 ) is uncertain for low- mass 
stars below about 3,700 K effective temperature. 



Based on the calc ulations by ISelsis et al.l (|2007l ). 
Mavor et al.l (120091 ) concluded that GL 581 d is a 
habitable planet, considering that the revised orbit 
means that it receives more than 30 % more stellar 
energy than previously thought. Because according 
to our study the low CO2 and most of the medium 
CO2 scenarios are uninhabitable (see Figs. [5] and ED, 
this disagrees with their general implication. This 
can be motivated in more detail as follows: Given 
that the equilibrium temperature of a planet only 
increases as the fourth root of stellar energy input, 
the increase of received stellar energy would only 
lead to an increase of about 10 K in equilibrium 
temperature (from about 185 K to 195 K). Thus, 
only an atmosphere providing a large greenhouse 
effect would be able to warm the planet above 273 
K, hence result in habitable conditions (see for ex- 
ample Fig. E]). 

The recent atm o spher ic modeling studies by 
Wordsworth et al. feoioh qualitatively a gree with 
our simulations. Wordsworth et al. ( 20ld ) use a ID 
radiative-convective model like ours. They also in- 
corporate a correlated-k approach for the radiative 
transfer, similar to what what we did. However, 
they used a different parameterization of the contin- 
uum absorption of CO2 than the one used here (for 
a sensitivity study of CO2 continuum absorption, 
again see Appendix [B]). Their calculated surface 
temperatures for the high-C02 cloud-free cases (set 
G3 in Table [3|) are comparable with the results ob- 
tained here (about 310 K for a 10 bar atmosphere, 
about 350 K for a 20 bar atmosphere). 



6. Conclusions 

Detailed model calculations of possible atmospheres 
for the low-mass extrasolar planet GL 581 d have 
been presented in this explorative modeling study. 
Using an improved ID radiative-convective climate 
model, several key atmospheric parameters (e.g. 
surface pressure, atmospheric composition) were 
varied to investigate their influence on resulting 
surface conditions. The planetary scenarios inves- 
tigated here are consistent with assumptions made 
in the literature regarding surface pressures and at- 
mospheric compositions of terrestrial planets. 
GL 581 d is a potentially habitable planet, because 
for massive CO2 atmospheres (5 or more bar sur- 
face pressure with CO2 concentrations of 95 %, 20 
bar with 5 % CO2), the surface temperatures ex- 
ceeded 273 K, i.e. the freezing point of water. For 
these massive CO2 atmospheres, the surface tem- 
peratures could be as high as 357 K. 
The results of our high CO2 scenarios are con- 
firmed by the very r ecent model results of 
Wordsworth et al. (|20ld ). But our model calcula- 
tions imply that also atmospheres that are not CO2- 
dominated (i.e., 5% vmr instead of 95% vmr) could 
result in habitable conditions for GL 581 d. 
For atmospheric scenarios with less CO2, the planet 
was found to be uninhabitable in our calculations. 
Nevertheless, GL 581 d is the first extrasolar 
(potentially terrestrial) planet where habitable 
conditions are at least conceivable within a 
reasonable range of surface pressures and CO2 
concentrations. 
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Appendix A: Comparative studies used for 
model validations 

A.l. Runs 

Several sets of comparison runs were performed to 
compare the model results of the improved model 
with other published work and benchmark calcula- 
tions. 

For all comparative runs, the central star was the 
Sun. The solar input spectrum is based on the high- 
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resol ution spectrum provid ed by lOuevmardl \2004 ) Table A.l. Runs performed to test the sensitivity 
(e.g., Kitzmann et al. 2010l ). of the model to variations of CO2 content 



The first set of comparison runs tested the response 
of the model to doubling/quadrupling the CO2 con- 
tent within the modern Earth atmosphere to com- 
pare the results with investigations regarding the 
anthropogenic greenhouse effect (runs W1-W8 in 
Table GO}. 



The second set of comparison runs that tested the 
model sensitivity to changes in CO2 content was 
performed with a reduced solar luminosity S of 0.8 
times the present value to compare to published 
calculations of different models with the same as- 
sumptions (runs S1-S6 in Table |A~T]) . 



Run 


S 


Pn-j [barl 


PCO2 


[barl 

L J 


Wl 


1 


0.77 


10" 4 




W2 


1 


0.77 


3.55 


• 10~ 4 


W3 


1 


0.77 


7.1 • 


10" 4 


W4 


1 


0.77 


10" 3 




W5 


1 


0.77 


1.42 


• 10~ 3 


W6 


1 


0.77 


2.84 


• 10" 3 


W7 


1 


0.77 


5.68 


• 10" 3 


W8 


1 


0.77 


10- 2 




SI 


0.8 


0.8 


10" 3 




S2 


0.8 


0.8 


10~ 2 




S3 


0.8 


0.8 


1.5 • 


10-2 


S4 


0.8 


0.8 


10- 1 




S5 


0.8 


0.8 


1. 




S6 


0.8 


0.8 


10.0 





The simulation to which the runs from Table 
I A. 1 1 were compared included an adjustment of 
the model surface albedo in a way that mod- 
ern Earth reference calculations yield a surface 
temperature of 288 K. For the model used here, 
however, a reference surface temperature of 284.5 
K was chosen to account for the missing green- 
house effect of methane, ozone, and nitrous ox- 
ide, which are not included in the present model. 
This value of the surface tempera ture is calcu 



lated with a modern Earth model (Grenfell et al 



2007bl . when using the RRTM radiative transfer 



scheme developed by iMlawer et al.l 119971 ) upon ex- 
cluding the greenhouse effect provided by methane, 
ozone, and nitrous oxide. The model surface albedo 
is 0.24, which is slightly higher than the value of 
0.21 from von Paris et al. ( 20081 ). but still compati- 
ble with surface albedos of other cloud-free mod- 



els in the literature (e .g., iGoldblatt et al.l l2009al . 



Haqq-Misra et al. 20081 . both studies using 0.23 



as their surface albedo value). Note that the ac- 
tual global value for Earth is approximately 0.13 
( Kitzmann et al. 2010l . Rossow <fe Schiffer 19991 ). as 
stated above. 



The third set of comparison runs followed a pro- 
posed evolutionar y sequ e nce o f the Earth's atmo- 
sphere based on iHarti (|l978l ) in terms of time 
if, before present (runs H1-H8 in Table IA.2p . 
To this end, the model surface albedo was ad- 
justed to yield a surface temperature of 288 K 
for present luminosity to be consistent with model 
studies using this sequence, which also introduced 
this particular adjustment feasting &: Ackerman 



19861 . iKiehl & Dickinsonl Il987h. The resulting sur- 
face albedo was 0.185, in contrast to the assumed 
0.24 for all other runs. The assumed N2 partial pres- 
sure was 770 mb. 
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Table A. 2. Runs for the evolution sequence of lHart 
(jl978h 



Run 


h [Gy] 


S 


c C o 2 [ vmr 


HI 


0.0 


1.0 


3.3 10" 


-4 


H2 


0.5 


0.972 


3.3 10" 


-4 


H3 


1.0 


0.944 


6.5 10" 


-4 


H4 


1.5 


0.917 


2.9 10" 


-3 


H5 


2.0 


0.889 


8.6 10" 


-3 


H6 


2.5 


0.861 


1.8 10" 


-2 


H7 


3.0 


0.833 


3.3 10" 


-2 


H8 


4.25 


0.764 


0.31 





20 



15 



10 



E 



5 - 



-_ 

0.01 



□ Goldblatt et al. 2009 
x This work, RRTM 
o This work, MRAC 



0.10 



.00 



[mb] 



I 0.00 



Fig. A.l. Response of surface temperature to in- 
creases in CO2: Temperature change with respect 
to run Wl . 

A. 2. Results from comparative studies 



Figure |A~T1 summarizes the results of the runs of set 
W (Table IA.1|) . i.e. tests for the assumed anthro- 
pogenic greenhouse effect. Shown is the increase 
AT of surface temperature with respect to run Wl, 
which is a measure of the strength of the greenhouse 
effect caused by CO2. 



Figure IA.1I also shows values for the same numer- 
ical experiment performed with two other mod- 
els. The first model is actually a version of the 



model used here, but it uses RRTM (Mlawer et al 



19971 ) instead of MRAC (see model description 
above) in the IR radiative trans f er. Th e sec- 
ond model is from iGoldblatt et all rt20091oh (also 



used in IGoldblatt et al.1 l2009al ). which is also a 



ID radiative-convective model. It differs from the 
model used here in that the correlated-k radia- 
tive transfer in the IR is based on Hitran 1992. 
Furthermore, the radiative transfer uses a 2-stream 
method instead of the diffusivity approximation 
used by MRAC. The main difference though is the 
numerical scheme employed to re a ch the steady 



state atmosphere. IGoldblatt et~aH (|2009bh 



use a 



Newton- Raphson method, not a time-stepping al- 
gorithm. Figure IA.1I shows that the agreement be- 
tween the three models is indeed very good. 



Figure IA.2I compares the calculations for set 
S from Table jA.l[ _with calculations done by 



Haqq-Misra et al. ( 20081 ) (similar to runs S3 and 



S4). They used an updated version (e.g., new ab- 
sorption coefficients) of the ID r adiati ve-convective 
model of Kasting Sz Ackerman ( 19861 ). The differ- 
ence between their model and ours is that the 
IR radiative transfer is performed by exponential- 
sum fitting of transmission instead of correlated- 
k. Figure I A. 21 implies that the results from 



Haqq-Misra et al. ( 20081 ) and this work for the two 



specific runs S3 and S4 agree very well. 
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Fig. A. 2. Response of surface temperature to in- 
creasing CO2 partial pressures at a reduced solar 
constant (0.8 times present-day value; set S from 
Table fATTl) 

The results of the model calculations for set 
H from Table IA.2I were compare d with the 
model results of iKiehl &z Dickinson ( 1987 ) and 
Kasting fc Ackermanl (119861 ) (see F igs. 7 and 8 in 



Kiehl fc Dickinson! 1 19871 : Fig. 2 in iKastind 119891 ) 



and are shown in Fig. IA.31 All models included in 
this test are ID radiative-convective models. 
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Fig. A. 3. Surface tempe rature s calcu lated for the 
evolutionary sequence of Hart ( 1978I ). set H from 
Table lA2l 



The model by Kasting fc Ackerman ( 1986I ) uses 
band models and exponential sums in the radia- 
tive transfer instead of the correlated-k. Also, they 



calculated the relative humidity in the atmosphere 
by a slightly different approach than the one used 
by us. 



The study by IKiehl fe Dickinson! ()1987l ) uses a rel- 
atively high-resolution radiative transfer code (500 
intervals in the thermal region) based on band mod- 
els. A second difference is the convective adjustment 
scheme that is not based on a simple formulation 
of the adiabatic lapse rate, the approach used by 
Kasting fc Ackermanl ( 19861 ) and here. 
Calculated surface temperatures are within 2 K of 
the results of the two other studies. These differ- 
ences are quite small given that the models differ 
with respect to e.g., radiative transfer and the treat- 
ment of convection, as stated above. 

A. 3. Conclusions from comparative studies 

The purpose of this section was to test the up- 
dated and improved model against other pub- 
lished model results. It was shown that the model 



agrees well with recent mod e ls of von Paris et al 



(l2008h and lHaaa-Misra et al.1 (hoosh as well a s ear 
lier model calculations dKiehl &: Dickinson 1987 : 



Kasting fc Ackermanl Il986l ). Furthermore, by com- 
paring with benchmark studies regarding the an- 
thropogenic greenhouse effect, the model was found 
to ag ree also very well wit h other model results 
(e.g., iGoldblatt etaUl2009bl ). 



Appendix B: Sensitivity studies for the 
effect of radiative transfer approximations 
on surface temperature 

In order to assess the sensitivity of the surface tem- 
peratures calculated for the atmospheric scenarios 
of Table [3] to details in the radiative transfer (such 
as continuum absorption or line mixing), an addi- 
tional sensitivity study was performed. Based on 
this sensitivity study, the results regarding surface 
temperatures and habitability can be interpreted 
with more confidence. Also, it serves as an illus- 
tration of the uncertainties associated with atmo- 
spheric modeling in general. 

The sensitivity runs have been performed for the 
simulations with medium and high CO2 concentra- 
tions (sets G2-G3 in Table |3|), because these runs 
resulted in habitable surface conditions. 
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B.l. Performed runs 

Both H2O and CO2 show significant collision- 
induced continuum absorption in the mid- to far- 
infrared. Especially the formulations of the foreign 
continua of H2O and CO2 could be relatively un- 
certain, because they are calcu lated for N2-O2 bac k- 
ground atmospheres (see, e.g. Halevv et al. 20091 ) . 
The first set of sensitivity runs was performed in or- 
der to assess the influence of the H2O continuum on 
the surface temperature (runs CH1-CH4 in Table 
IB.lj) . These tests were done for the 5, 10, and 20 
bar atmosphere of set G3 as well as the 20 bar run 
from set G2 (Table [3]) by multiplying the H2O con- 
tinuum by arbitrary factors /h 2 o of 0.2 and 0.5. 
The second set of sensitivity runs was performed to 
assess the influence of the CO2 self-continuum on 
the surface temperature (runs CS1-CS4 in Table 
IB.lj) . The CO2 self-continuum was multiplied by 
arbitrary factors /co 2 of 0.2 and 0.5, respectively. 
To test the influence of CO2 foreign-continuum ab- 
sorption, the continuum was removed for the same 
runs, i.e. multiplied by (runs CF1-CF4 in Table 
IB.lj) . These tests were performed for the same runs 
as the H2O sensitivity runs. 



An additional challenge in the radiative trans- 
fer arises from so-called line-mixing. Assuming a 
Voigt line profile (i.e. a convolution of Lorentz and 
Gauss profiles) is no longer justified. Comparisons 
of computer sim u lation s with experimental data by 
Rodrigues et al. (jl999h have shown that including 



line-mixing into the calculations can result in a de- 
crease of absorption coefficients by up to a factor of 
2. Line-mixing paramet ers for CO2 are includ ed in 
the HITRAN database dRothman et all 12005 1. but 
not in the HiTemp database (jRothman et al.lll995l ) 
used for the absorption coefficients in the model. 
Hence, a Voigt profile was assumed for all lines at 
all pressures dur ing the calculations o f the absorp- 
tion coefficients ( von Paris et al. 20081 ). In order to 
test the sensitivity of the results to line-mixing, a 
third set of sensitivity runs was performed (runs 
LM1-LM4 in Table IBTT]) . Here, the optical depth in 
the main atmospheric band of CO2 (15 /im) was 
arbitrarily decreased in the lower troposphere (i.e. 
at pressures higher than 100 mb) by factors /lm 
of 0.2 and 0.5, respectively. Again, these tests were 
performed for the same runs as the H2O sensitivity 
runs. 

Table IBTT1 summarizes the sensitivity runs. 



Table B.l. Sensitivity runs performed for GL 581 
d (/h 2 H2O continuum; fco 2 CO2 continuum; /lm 
line mixing factor) 



Set 


Control Run 


/h 2 o 


/co 2 


/lm 


CHI 


G3 20 bar 


0.2,0.5 


1 


1 


CH2 


G3 10 bar 


0.2,0.5 


1 


1 


CH3 


G3 5 bar 


0.2,0.5 


1 


1 


CH4 


G2 20 bar 


0.2,0.5 


1 


1 


CS1 


G3 20 bar 


1 


0.2,0.5 


1 


CS2 


G3 10 bar 


1 


0.2,0.5 


1 


CS3 


G3 5 bar 


1 


0.2,0.5 


1 


CS4 


G2 20 bar 


1 


0.2,0.5 


1 


CF1 


G3 20 bar 


1 





1 


Cf2 


G3 10 bar 


1 





1 


CF3 


G3 5 bar 


1 





1 


CF4 


G2 20 bar 


1 





1 


LM1 


G3 20 bar 


1 


1 


0.2,0.5 


LM2 


G3 10 bar 


1 


1 


0.2,0.5 


LM3 


G3 5 bar 


1 


1 


0.2,0.5 


LM4 


G2 20 bar 


1 


1 


0.2,0.5 



B.2. Effect of H2O continuum 

The runs of sets CH1-CH4 from Table IB. II were 
done in order to test the influence of the H2O con- 
tinuum on calculated surface temperatures. The ob- 
tained results did not differ significantly from the 
control runs. For a decrease of H2O continuum of 
a factor of 5, the corresponding decrease of surface 
temperature was less than 2 K in each case. 
Thus, these results indicate that the formulation of 
the H2O continuum is not critical for the calcula- 
tions presented here. 

B.3. Effect of CO2 continuum 

Sensitivity tests were performed to assess the influ- 
ence of the CO2 continuum on surface temperature 
(sets CS1-CS4 and CF1-CF4 from Table EB). For 
the CO2 self-continuum, the effect was relatively 
large. The strong mid-IR band around 7 /im dom- 
inates absorption and is the main opacity source 
in this spectral region. Also, the far-IR bands show 
significant contribution to the overall absorption. 
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Hence, the effect of a decrease of the CO2 contin- 
uum can be expected to be very strong. 
For set CS1, the surface temperature decreased 
from 357 K (control run) to 346 K and 329 K on 
decreasing the continuum absorption by a factor of 
2 and 5, respectively. 

Table IB.2I summarizes the effect on surface temper- 
atures. 



Table B.2. Surface temperatures (in K) for the 
sets CS1-CS4 from Table EE 



Set Control Run /co 2 =0-2 /co 2 = 0-5 



CS1 
CS2 
CS3 
CS4 



357 
322 
287 
313 



329 
292 
257 
282 



346 
311 
274 
300 



From the results in Table IB.2| one infers that the 
uncertainties in the CO2 self continuum opacity will 
not alter the principle conclusions on surface hab- 
itability provided that they do not exceed a fac- 
tor of 2-5. The results of the sensitivity tests im- 
ply, however, that more detailed modeling and mea- 
surements of the CO2 self-continuum absorption are 
needed in the future. 

Upon removing the CO2 foreign continuum (sets 
CF1-CF4 in Table iBTTj) . surface temperatures de- 
crease by 1-3 K for the high C0 2 runs (CF1-CF3). 
For the 20 bar medium CO2 run, however, exclud- 
ing the CO2 foreign continuum decreased the sur- 
face temperature by 33 K to 280 K. This implies 
that the foreign continuum is an important opacity 
source and should be included in all future simula- 
tions. 

B.4. Effect of line mixing 

The effect of line mixing on surface temperatures 
was investigated with the sets LM1-LM4 from Table 

EU 

For all cases, the surface temperatures were al- 
most unaffected (less than 1 K decrease). These ef- 
fects are important for interpreting measurements 
on Mars or Venus, but are not likely to significantly 



affect results regarding exoplanets, where only first- 
order estimates can be done so far. 

B.5. Summary of sensitivity studies 

The purpose of this section was to investigate the 
influence of the radiative transfer approximations 
on the calculated surface temperatures. 
It was shown that the model results do depend to a 
certain extent on the details of the radiative trans- 
fer. However, only the effect of the CO2 continuum 
is likely to be important for our study. 
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